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Abstract. We give an efficient complete parametrization of wavelet matrices 
of rank 771, genus g + 1, and degree g, which are naturally identified with corre- 
sponding polynomial paraunitary matrix- functions. The parametrization depends 
on Wiener-Hopf factorization of unitary matrix-functions with constant determi- 
nant given in the unit circle. This method allows us to construct in real time the 
coefficients of wavelet matrices from the above class. 



A wavelet matrix A = 
infinite vectors 



(cfi) of rank m consists of m rows of possibly 
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(£ E C, satisfying following two conditions. 



(i) Quadratic condition: 
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fii) Linear condition: 



(3) 



E 



a j 



m5 



r.O 



where 5 stands for the Kronecker symbol. 

In this paper we assume that (1) is compact, i.e. only finite number of 
its entries are different from 0. Therefore, the series in (2) and (3) are only 
formally infinite, and no problem of convergence appears. 

The quadratic condition (2) asserts that the rows of a wavelet matrix 
a r := (Qj)j= _oo,oo have length equal to \pm and that they are pairwise 
orthogonal when shifted by an arbitrary multiple of m. The first row a is 
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called the scaling vector or low-pass filter, while remaining rows a r are called 
the wavelet vectors or high-pass filters. In signal processing applications, the 
linear constraint (3) implies that a constant signal emerges from the first 
subband of the maltirate filter bank (1). 

Associate to each wavelet matrix A the matrix function A(z) as follows: 
let Ak, k G Z, be submatrices of A of size mxm defined by A k = (a r km+s ), 
< s,r <m — 1, in other words, (1) is expressed in terms of block matrices 
in the form 

A=(--- ,A_ 1 ,A ,A 1 ,A 2 ,---), 

and assume 

oo 

(4) A(z) = £ A k z k . 

k=—oo 

Obviously, there is one-to-one correspondence between matrices (1) and 
formal series expansions (4), and, for a compact matrix, the corresponding 
matrix-function is a Laurent polynomial. 

It can be verified that the quadratic and the linear constraints on A are 
equivalent, respectively, to the following two conditions on A(z): 

(5) A(z)A*(z~ 1 ) = ml, 

where A*(z' 1 ) := YlkL-oo ^l z ^ k * s the adjoint of A(z), and 

m 

(6) y^ Ajj(l) = m^i, 1 < i < m, 

3=1 

where A(z) = (Ay(^)). ._ . The condition (5) means that A is a parauni- 
tary matrix-function. 

If U is a unitary matrix of size mxm, U G U(m), and A(z) satisfies (5), 
then UA(z) satisfies (5) as well. Furthermore, for each paraunitary matrix- 
function A(z), there exists and one can explicitly construct a unitary matrix 
U such that UA(z) satisfies the linear condition (6) as well. If U and U' are 
two such matrices, then 

"•-(Sty* 

where V G U(m—1). Thus the construction of paraunitary matrix-functions 
are decisive for construction of wavelet matrices. 

It is said that a wavelet matrix (1) has the rank m and the genus g, 
A G WM(m, g; C), if the corresponding matrix-function A(z) has a form 

(7) A(z) = J2A k z k . 

k=0 



It can be easily shown (see [5, p. 58] that the determinant of a paraunitary 
matrix-function A(z) is a monomial in z, that is, there is a nonnegative 
integer d, called the degree of A.(z), such that 

det A(z) = cz d . 

Generically, (7) has degree g — 1, although in specific degenerated cases, it 
can be larger or smaller than g — 1 . 

The relation between compact wavelet matrices and compactly supported 
wavelet systems as orthonormal functions in L 2 (R) is well-known (see [5], 
Ch. 5). 

Theorem ([1], [4], for rank 2; [5, pp. 87, 91], for rank m > 2): Let 

AeWM(m,g;C) 
be a wavelet matrix and consider the functional difference equation 

rag— 1 

(8) 4>{x) = 2_. (^^{mx — k) 

k=0 

called the scaling equation associated with A. Then, there exists a unique 
(j) E L 2 (R), called the scaling function, which solves (8) and satisfies 



(x) dx = 1 and supp <fi C 



0,(5-1) (-^-)+l 

m — 1 



Furthermore, if we define wavelet functions (associated with A) by the 
formula 

mg— 1 

ip r (x) = 2_. Ci r k ( t ) {' mx ~ &)) 1 < r < m, 

fc=0 

and consider the collection of functions 

4>jk(x) = m? ' 2 <\>{m? x — k), j,k E Z, 
V^ fc (£) = mP' 2 ilf{mPx — k), 1 < r < m; j,k E Z, 

called the wavelet system W[*4] (associated with wavelet matrix ^4), then 
there exists an /^-convergent expansion for each f E L2. 

00 m— 1 00 00 

(9) f(x)= Yl c kMx) + ^2Yl Yl cr jk^ r jk(x), 

k=— 00 r=l j'=0 k=— 00 



where the coefficients are given by 



Cfc = / f(x)(j)ok{x)dx, 



DO 
OO 



c r jk= / f(x)ifj r jk (x)dx. 

J —oo 

Remark For most wavelet matrices A, the wavelet system W[^4] is a 
complete orthonormal system and hence an orthonormal bases for L 2 (R), 
which would imply the above theorem. However, for some wavelet matrices, 
the system W[*4] is not orthonormal, and yet (9) is always true, which 
means that W[^4] is a tight frame. 

Independently from the above mentioned connection between the wavelet 
matrices and associated wavelet systems, the former can be directly used 
in various discrete signal processing applications. Namely, the following 
theorem is one of the key links between the mathematical theory of wavelets 
and its practical applications. 

Theorem (wavelet matrix expansion, [5, p. 80]): Let 

/: Z^C 

be an arbitrary function (discrete signal) and let 

A=(a r k )eWM(m,g;C) 

be a wavelet matrix of rank m and genus g. Then / has a unique wavelet 
matrix expansion 

m— 1 oo 
J\ n ) = /_^i Z-^i C k a mk+m 

where 
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C k ~ ~ /. J\ n ) a mk+n- 

m *-^ 
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The wavelet matrix expansion is locally finite; that is, for given n, only 
finitely many terms of the series are different from 0. 

From whatever said above it is evident the theoretical and practical impor- 
tance of deeper understanding an internal structure of paraunitary matrix- 
functions which would allow to construct efficiently a wide class of such 
matrices. So far, the only way of classification of paraunitary matrix func- 
tions was via the following factorization theorem. This theorem resembles 
the factorization of polynomials of degree d according to their d roots and 
highest coefficients. 

For a unit column vector v G C m , v*v — 1, let 

(10) V(z) := I-vv* + vv*z. 



Obviously, (10) is a polynomial matrix function of order 1. It can be shown 
that V(z) is the paraunitary matrix-function of degree 1 (see [5, p. 59]) and 
it is called primitive. 

Theorem (Paraunitary Matrix Factorization, [5, p. 60]): A parau- 
nitary matrix-function (7) of degree d, where A g ^\ ^ can be factorized 



as 



A(z) = V 1 (z)V 2 (z)---V d (z)U 

where Vj(z), j = 1,2, ... ,d, are primitive paraunitary matrix-functions and 
U is a (constant) unitary matrix. 

We propose absolutely new way of parametrization of paraunitary matrix- 
functions of rank m, genus g + 1 and order g, which depends on Wiener- 
Hopf factorization of unitary matrix-functions (with constant determinant) 
given on the unit circle in the complex plane. Actually this method was 
developed in [3] , [2] and it allows to construct efficiently matrix-functions of 
the above type (consequently, to prepare the coefficients of the whole class 
of compactly supported wavelets) in real time. 

Let 



A(z) 



/A lfl (*) A h2 (z) 
M,i( z ) A 2 ,2(*) 
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A(z)A*(z- 1 ) = I, 



det A(z) = cz 9 , 



\c\ = 1. 



We first convert the matrix-function A(z) into a unitary (on the unit circle) 
matrix-function U(z) by dividing any row of A(z) (say, the last row, to be 
specific) by z 9 : 
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Then we have 



U{z)U*{z)=I, for |3| = 1, 
det U(z) = c, \c\ = 1, 

1 < r < m; 1 < j < m, 
U m ,j e L g , 1 < j < m, 



U rj G L + g 



where 



L+ = {f:f(z) = J2c k z k } : 



k=0 
9 



L- = {f:f(z) = J2ckZ~ k }. 



fc=0 

Thus, the two theorems below give simple and transparent way of one- 
to-one parametrization of paraunitary matrix-functions of rank m, genus 
g + 1 and degree g. To compare this with the above presented factorization 
theorem from the simplicity point of view, the proposed parametrization 
resembles the classification of polynomials of degree d according to their 
d + 1 coefficients. 

Theorem 1. (see [2], p. 22) For each mxm matrix-function F(z) of form 

( 1 • • • 0\ 
1 ••• 
1 ••• 
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where 
(12) 
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¥j e L g , 1 < j < m- 1, 

there exists a unitary matrix-function U(z) (unique up to a constant unitary 
right multiplier) of form 

( Wll Wl2 • • • Ui m \ 

U 2 \ U 22 ■■■ U 2m 



(13) 

where 
(14) 
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with constant determinant, such that 

(15) F(z)U(z) E L+. 

Theorem 2. (see [3]) For each unitary matrix-function U(z), with constant 
determinant, of form (13), (14), there exists an unique F(z) of form (11), 
(12) such that (15) holds. 

Observe that if we denote by F_(z) the matrix function of type (11) where 
each ifij is replaced by — ipj, j — 1, 2, . . . , m — 1, then F_{z) = (F(z))^ 1 , so 
that the equation 

U(z) = F_{z) ■ F(z)U{z) 
gives the rigt Wiener-Hopf factorization of U(z). 

In the end, we should mention that less than 1 sc computer time is re- 
quired to compute coefficients of matrix-function (13) in Theorem 1 when- 
ever coefficients of functions ipj, j = 1, 2, . . . , m — 1, are selected in (11) for 
such large dimensions as m = 30 and g — 50. This speed of calculations 
opens possibility to choose the optimal wavelet matrix for specific problem, 
which is the most important step in practical applications, by total selection. 
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